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Using an exact numerical solution to Maxwell's equations, we 
determine the optimum refractive index profile parameter a for a 
germania- doped silica lightguide. Our results agree closely with the 
earlier work of Olshansky and Keck and with the work of Marcatili. 
However, differences exist that may be important in the manufacture 
of very high bandwidth lightguides. 

I. INTRODUCTION 

It is now well established that a lightguide with a near parabolic 
index profile can have single frequency bandwidths of 1 GHz -km or 
greater. 1-3 This is expected to happen if the profile is smooth, the 
central dip is narrow or small, the profile parameter is properly chosen, 
and modes that have a significant amount of their energy transported 
in the cladding are eliminated. 4 " 6 Our understanding of these factors 
has largely come from approximate solutions to Maxwell's equations 
for lightguides. The availability of computer software for exact numer- 
ical solutions 7 to Maxwell's equations allows us to reexamine some of 
these effects and to look at other factors that may influence lightguide 
bandwidths. In this paper, we determine optimum power law profile 
parameters (i.e., a parameters) as a function of wavelength for a 
germania-doped silica lightguide. 

II. PROCEDURE 

We study a lightguide with a 13.5 mole percent germania-doped 
silica core and a pure silica cladding. The core is assumed to be 
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perfectly cylindrical with a diameter of 50 microns. The cladding is 
considered to be essentially infinite. The doping profile is such that 
the refractive index is a power law function of distance R. Thus, 



N = Ni + (N 2 - Ni) 



ReeJ ' 



(1) 



where Ni is the refractive index at the center of the core, N 2 is the 
refractive index in the cladding, R cc is the radius of the core, and a is 
the profile parameter. We assume that this functional form holds at all 
wavelengths of interest. 

We have accurately determined the wavelength variation of the 
refractive index of silica and some doped silicas. Figure 1 gives the 
data for 13.5 mole percent germania-doped silica and pure silica. As is 
well known, the refractive index falls off slowly with wavelength for 
both of these glasses. 8 We express the data by the following modified 
Sellmeier formula: 9 

Nt = Co + CiA 2 + C 2 \ 4 + C 3 /(X 2 - 0.035) 

+ C 4 /(\ 2 - 0.035) 2 + C 5 /(X 2 - 0.035) 3 , (2) 



1.47 - 
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Fig. 1 — The wavelength dependence of 13.5 mole percent germania-doped silica and 
pure silica. 
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where A is in microns and i = 1, 2. Table I lists the coefficients for eq. 
(2) for both 13.5 mole percent germania-doped silica and silica. 

Our numerical procedure to solve Maxwell's equations has been 
described in detail in an earlier paper. 7 Thus, we only give a very brief 
description here. The method is similar to that described by Vigants 
and Schlesinger, 10 and Vassel. 11 We write Maxwell's equations in 
cylindrical coordinates as 

^ = p- l A( P )T(p), (3) 

up 

where T{p) is a column vector whose four elements are related to the 
tangential components of the electromagnetic field. They are as fol- 
lows: 

r, = e 2 , 

T 2 = -ipH+Zo, 

(4) 

r 3 = pE+, 

T4 = —ffizZo, 

where Zo is the wave impedance of free space and the variable p is 
defined as KoR. 

The 4x4 matrix A(p) can be written as 



(5) 



This matrix describes the properties of the media. 

We solve eq. (3) by an optimized fourth-order Runge-Kutta proce- 
dure. 12 This is a higher-accuracy formula than used in our earlier work. 
The most important computational results are the effective indices N e 
for the various modes in the lightguide. 

The group indices N M are calculated from the effective indices by 
the following formula: 

N e = N e - XdN e /dX. (6) 

Table I — Coefficients for eq. (2) 






m/K) - 1 





-MNe/K 


P 2 k-M 2 





MN e 








MNe/K 





p 2 - (M 2 /k) 


-MN e 





N 2 e - K 






13.5 Mole Percent Germania- 

Doped Core Silica Cladding 

Co = 1.4706868 1.4508554 

C, = -0.0026870 -0.0031268 

C 2 = -0.0000356 -0.0000381 

C 3 = 0.0035756 0.0030270 

C 4 = -0.0000828 -0.0000779 

C 5 = 0.0000018 0.0000018 
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The derivative is obtained from numerical calculations of N e at three 
closely spaced wavelengths. 

For a lightguide with the materials and profile described earlier, a 
simple approximate relationship exists between group index N g and 
effective index N e for certain modes. In particular, for a given profile 
parameter a, angular mode number M, wavelength A, and for modes 
far from cutoff, the effective indices are nearly linearly related to the 
group indices. Evidence for this can be obtained from Fig. 2, curve (a). 
In this case, we have X = 0.6328 micron, M = 0, and a = 2. The solid 
line is the linear least-squares fit to the data when the last two modes 
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Fig. 2 — Plot of group index (2V») versus effective index {N e ) at a wavelength (X) of 
0.6328 micron with M = 0. In (a) the profile parameter a is 2.0 while in (b) o is 2.2. Note 
that a good linear relationship exists between N g and N e if the last few modes are 
neglected. 



458 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1981 







CORE INDEX 




CLADDING 
INDEX 


1.4900 






2.40 
2.30 


T 

1 ns/km 
i 


1.4895 








2.20 


1 .4890 








-2.15 

-2.10 
- 2.08 


1 .4885 








-2.05 
-2.00 


1 .4880 










1.4875 


- 


X = 0.82 




-1,90 
-1.80 


1.4870 


- 








1.4865 


1 


1 


1 


1 



1.465 



1.455 



Fig. 3 — Plot of group index versus effective index for X = 0.82 micron and M = 0. 
When a = 2.08 the group indices are all nearly the same except for those modes being 
influenced by the cladding. 



are deleted. (The last mode is too far off scale to be plotted.) Obviously 
a nearly linear relationship exists between N e and N e . In addition, 
there is a very slight systematic deviation from the least-squares line. 
This suggests that additional terms are needed to completely describe 
the data. However, this will be of no consequence to the analysis to be 
described here. 

The last two modes (those that deviate strongly from the line) have 
a significant fraction of their energy being transported in the cladding. 
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Thus, they are not tightly bound to the core and normally are not 
considered in an optimum a calculation. 

An optimum profile parameter a, for a particular family of modes, 
is one that minimizes the spread in group indices for that family at a 
particular wavelength. For example, we might consider all the modes 
associated with a specified M value. In this case, the optimum a would 
be the one that gives the smallest slope to the least-squares line 
relating N e and N g . 
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Pig. 4 — Part of the data of Fig. 3 on an expanded scale. 
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Fig. 5— Plot of group index versus effective index for X = 1.06 microns and M = 0. 

Figure 2, curve (b) shows a near optimum value of a = 2.2 for M = 
and X = 0.6328 micron. We see that the least-squares line is now 
nearly horizontal and that evidence exists for the necessity of higher- 
order terms to exactly fit the data. The strong drop-off of the last few 
modes is still evident. 

A very good overall optimum a for all M values will be nearly 
identical to the optimum a for M = 0. As will be shown shortly, this 
results in excellent equalization of the lower-order modes and a good 
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equalization of the higher-order modes. For most engineering appli- 
cations, optimum a's in the wavelength range of 0.63 to 1.55 microns 
are needed. 

III. RESULTS 

Figures 3 through 7 show plots of N g versus N e for various a 
parameters and wavelengths for M = 0. Figure 4 shows part of the 
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Fig. 6 — Plot of group index versus effective index for X = 1.32 microns and M = 0. 
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Fig. 7— Plot of group index versus effective index for X = 1.55 microns and M = 0. 

data of Fig. 3 on an expanded scale. Note that in all a nearly linear 
relationship exists between N e and N e , except for the last couple of 
modes. As mentioned before, these modes are nearing cutoff and are 
usually not considered in an optimum a determination. 

At each wavelength, we see that there is an a value, a opt , that makes 
the slope of the line nearly zero. Table II lists «o P t values for wave- 
lengths of engineering importance and also for a few other wavelengths 
as well. 

The data of Figs. 3-7 pertaining to the infrared wavelengths are 
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Table II — Optimum alpha values 



A (microns) 


Qopt 


0.6328 


2.211 


0.73 


2.130 


0.82 


2.081 


0.90 
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2.000 
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1.32 


1.884 


1.40 


1.856 


1.55 


1.807 



summarized in Figs. 8 and 9. Figure 8 shows the slope of the line 
relating N g to N e as a function of wavelength for various a values. We 
see that most infrared wavelengths of interest require an a value in 
the range from 1.8 to 2.1. Figure 9 shows the same data, but this time 
as a function of a for various wavelengths. The same conclusion is 
obvious. 

Figure 10 displays a calculation of the group indices for all the bound 
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Fig. 8 — Plot of the slope of N g versus TV, as a function of A for various a values. 
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modes in the fiber for X = 0.82 micron and a op t = 2.081. This plot 
gives a vivid representation of the fiber characteristics. We see that 
indeed a = 2.081 gives a good equalization for all the modes if those 
being influenced by the cladding are ignored. A delay-time calibration 
corresponding to 1 ns/km is also included in the plot. 



IV. DISCUSSION 

Within the framework of the WKB approximation, Marcatili 13 has 
obtained an exact result for a op t that minimizes pulse widths. His 
formula is as follows: 



where 



«o P t=(2-p)(l+ Vl-2A)-2, 



P = 



N>_ X rfA 
N° g ' A ' d\ ' 



(7) 



(8) 
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Fig. 9 — Plot of the slope of N M versus N e as a function of a for various A values. 
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Fig. 10 — Plot of group index versus effective index for all the bound modes in the 
lightguide. We see that for A = 0.82 micron and a = 2.081 there is a good equalization of 
all the group indices. 



Nl is the group index on axis, and A = 



N1-N2 
Ni 



(9) 



Also, within the same WKB framework, Olshansky and Keck 14 have 
obtained an approximate formula for a opt that rninimizes the rms width 
of a pulse. Approximations that coincide with Olshansky and Keek's 
can also be obtained from the exact result of Marcatili. 

It is clear that a op t depends upon the type and number of approxi- 
mations employed, how the various modes are weighted, and whether 
we wish to minimize pulse width, rms width, or some other parameter. 
Another possibility would be to define an optimum a that would 
include the effects of those modes having significant energy trans- 
ported in the cladding. This is likely to lead to an a quite different 
from that obtained in this work. 
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Figure 11 shows the results of Olshansky, Marcatili, and ourselves. 
There is great similarity in the shape of the curves and in the actual 
numbers. The deviation between our results and the WKB results is 
largest at the longer wavelengths. It is obvious that the agreement 
between all results is good and, at the present state of the manufac- 
turing art, can be considered the same. However, this may not neces- 
sarily be the case in the future. To illustrate (see Fig. 12), we calculate 
N g as a function of N e at 1.55 microns for M = and a opi = 1.775 
(Olshansky) and a op t = 1.807 (this work). There is an obvious difference 
in slope; this difference would be significant in the design of high 
bandwidth lightguides. 

A good representation of our optimum a values is given by the 
following formulas: 



0.63 micron a£ X s* 1.00 micron, 
aopt = 2.9970857 - 1.6647237A -I- 0.6629181A 2 , 



(10) 
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Fig. 11 — Plots of optimum profile parameter a versus wavelength for a 13.5 mole 
percent germania-doped silica lightguide. The index dispersion data employed is from 
Fleming's work. A clear similarity exists among the results of Marcatili, Olshansky, and 
this work. 
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Fig. 12 — A plot of N„ versus N e for A = 1.55 microns. a op t = 1.775 is calculated from 
the formula of Olshansky et al., and o op t = 1.807 is from this work. 



1.00 micron «s A s£ 1.55 microns, 
aopt = 2.3835734 - 0.4170708A + 0.0290252A 2 . 



(11) 



It is reasonable to assume that by linear interpolation or extrapola- 
tion we can generate the refractive index dispersion curves for ger- 
mania- doped silica in the range from to 16 mole percent. Using this 
data, we can then calculate the dependence of N g on N e for a number 
of levels of doping. Figure 13 shows a family of curves for M = 0, a = 
1.884, and A = 1.32 microns. There is a good equalization of N g for all 
the compositions. We also see that, as expected, the number of modes 
decreases as the doping level drops. Thus, at the 1 percent doping level 
for M = only one N e exists. As usual, the last couple of modes for 
each doping level show a substantial cladding effect. 

Finally, the data displayed in Figs. 8 and 9 may have some implica- 
tions to lightguide engineering. From Fig. 8, we see that the various a 
curves pass through zero with about the same slope. This suggests that 
the range of operating wavelengths will be the same for all design 
wavelengths. Likewise from Fig. 9, we see that the various X curves 
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Fig. 13 — A family of plots of N g versus N e for various doping percentages of germania. 
The a value is 1.884 and M = 0. 



pass through zero with similar slopes. This suggests that the allowable 
error in a will be about the same for all design as. 

V. ACKNOWLEDGMENTS 

We greatly appreciate the enthusiasm for this work expressed by M. 
I. Cohen, L. S. Watkins, and R. J. Klaiber. 



NUMERICAL CALCULATION OF OPTIMUM a 469 



REFERENCES 

1. D. B. Keck and R. Bouillie, "Measurements on High-Bandwidth Optical Wave- 

guides," Opt. Commun., 25 (April 1978), pp. 43-8. 

2. M. Nakahara et al., "Ultra Wide Bandwidth V.A.D. Fibre," Electron. Lett., 16 (May 

1980), pp. 391-2. 

3. M. Nakahara et al., "Fabrication of Low-Loss and Wide Bandwidth V.A.D. Optical 

Fibres at 1.3 urn Wavelength," Electron. Lett., 16 (January 1980), pp. 102-3. 

4. D. Marcuse, "Calculation of Bandwidth From Index Profiles of Optical Fibers 1: 

Theory," Appl. Opt., 18 (June 1979), pp. 2073-80. 

5. R. Olshansky, "Propagation in Glass Optical Waveguides," Rev. Mod. Phys., 51 

(April 1979), pp. 341-67. 

6. S. Ishikawa, K. Furuya, and Y. Suematsu, "Vector Wave Analysis of Broadband 

Multimode Optical Fibers With Optimum Refractive Index Distribution," J. Opt. 
Soc. Am., 68 (May 1978), pp. 577-83. 

7. G. E. Peterson et al., "An Exact Numerical Solution to Maxwell's Equations for 

Lightguides," B.S.T.J., 59 (September 1980), pp. 1175-96. 

8. J. W. Fleming, "Material Dispersion in Lightguide Glasses," Electron. Lett., 14 (May 

1978), pp. 326-8. 

9. M. Herzberger, Handbook of Physics, New York: McGraw-Hill, 1958, E. U. Condon 

and H. Odishaw, Eds., pp. 6-42 and 6-43. 

10. A. Vigants and S. P. Schlesinger, "Surface Waves on Radially Inhomogeneous 

Cylinders," IEEE Trans., MTT-10 (September 1962), pp. 375-82. 

11. M. 0. Vassel, "Calculation of Propagating Modes in a Graded Index Optical Fiber," 

Opto-Electron., 6 (July 1974), pp. 271-386. 

12. Anthony Ralston, "Runge-Kutta Methods With Minimum Error Bounds," Math. 

Comput., 16 (October 1962), pp. 431-7. 

13. E. A. J. Marcatili, "Modal Dispersion in Optical Fibers With Arbitrary Numerical 

Aperture and Profile Dispersion," B.S.T.J., 56, No. 1 (January 1977), pp. 49-63. 

14. R. Olshansky and D. B. Keck, "Pulse Broadening in Graded-Index Optical Fibers," 

Appl. Opt., 15 (February 1976), pp. 483-91. 



470 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1 981 



